Hyperfine tensors of nitrogen-vacancy center in diamond from ab initio calculations 
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We determine and analyze the charge and spin density distributions of nitrogen-vacancy (N-V) center in di- 
amond for both the ground and excited states by ab initio supercell calculations. We show that the hyperfine 
tensor of N nuclear spin is negative and strongly anisotropic in the excited state, in contrast to previous models 
used extensively to explain electron spin resonance measurements. In addition, we detect a significant redistri- 
bution of the spin density due to excitation that has serious implications for the quantum register applications of 
N-V center. 



Nitrogen-vacancy (N-V) centers in diamond have numer- 
ous peculiar properties that make them a very attractive solid 
state system for fundamental investigations of spin based phe- 
nomena. Recently, this defect has been proposed for several 
applications, like quantum information processing ||jl|2l|3|,|4J], 
ultrasensitive magnetometer lUllil], and measurement of zero- 
point fluctuations or preparing quantum-correlated spin-states 
over macroscopic distance [7]. In these measurements, a room 
temperature read-out of single nuclear spins in diamond has 
been achieved by coherently mapping nuclear spin states onto 
the electron spin of a single N-V center yl |8fl which can 
be optically polarized and read out with long coherence time 
JHOl. Particularly, this has been the basis in realization of a 
nuclear-spin-based quantum register 111 ill and multipartite en- 
tanglement among single spins at room temperature fczll . The 
polarization of a single nuclear spin has been achieved by us- 
ing either a combination of selective microwave excitation and 
controlled Larmor precession of the nuclear-spin state 111 ill or 
a level anticrossing in the excited state lfl3ll . Understanding 
the spin states and levels is of critical importance for optical 
control of N-V centers in both the ground and excited states. 
Especially, the hyperfine interaction couples the electron spin 
and nuclear spin, thus determination of hyperfine tensors of 
the nuclei with non-zero nuclear spin plays a key role both in 
creation of entaglement states and in the decoherence process 



Recently, the hyperfine signals in the ground IU3L Il5ll and 
excited [ 16] states have been detected in N-V centers but with 
contradicting interpretations. In a conventional electron para- 
magnetic resonance (EPR) spectrum on the ensemble of N- 
V centers, the 15 N signal was assumed to be positive with 
slight anisotropy in the ground state while Fuchs et al. and 
Jacques et al. assumed an isotropic negative hyperfine con- 
stant for 15 N in the ground state IU3[ llql based on previ- 



ous EPR and optically detected magnetic resonance (ODMR) 
measurements 11171 fl8ll . Recently, Fuchs et al. have reported 
that the hyperfine splitting of 15 N should be ^20 x larger in 
the excited state than in the ground state Jlrjll . The excited- 
state 15 N hyperfine signal was assumed to be isotropic in their 
model [16]. While we already addressed the hyperfine tensor 



of N in the ground state 



ady 

m 



the lack of the detailed study 



on & N hyperfine signal and the proximate 3 C isotopes both 
in the ground and excited states prohibits the understanding of 



the intriguing physical properties of this defect. In this Letter, 
we thoroughly investigate the hyperfine tensors of proximate 
15 N and 13 C isotopes of the N-V center both in the ground 
and excited states by means of high level ab initio supercell 
all-electron plane wave calculations. In addition, we analyse 
the overall charge and spin density distributions before and 
after the optical excitation. We show that the hyperfine con- 
stants of 15 N is positive and possibly slightly anisotropic in 
the ground state while negative and strongly anisotropic in 
the excited state. In addition, the hyperfine splittings of the 
proximate 13 C isotopes change also significantly that has seri- 
ous implications both in the interpretation of the excited-state 
spectroscopy signals and in the quantum-information applica- 
tions. 

The negatively charged NV-center in diamond lEoil consists 
of a substitutional nitrogen atom associated with a vacancy at 
an adjacent lattice site (Fig. QJ). The ground state has 3 ' Ai 
symmetry where one a\ defect level in the gap is fully occu- 
pied by electrons while the double degenerate e-level above 
that is occupied by only two electrons with parallel aligment 
of spins (Fig.[TJ)). Thus, this defect has S=l high spin ground 
state 02 ill . By promotion one electron from the ai-level to the 
e-level will result in the excited 3 E state. Both states can be 
described by conventional density functional methods. The 
ground state can be readily described by spin density func- 
tional theory while the excited state can be obtained by con- 
strained occupation of states 01911 . We optimized the geome- 
try for both the ground and excited states, and calculated the 
charge and spin densities at the optimized geometries. 

We applied the PBE functional B24ll to calculate the spin 
density of the defect. First, the diamond primitive lattice was 
optimized, then a simple cubic 512-atom supercell was con- 
structed from that. Finally, we placed the negatively charged 
NV-defect into the supercell, and optimized under the given 
electronic configuration. We utilized the VASP code for ge- 
ometry optimization |25]. We applied plane wave basis set 
(cut-off: 420 eV) with PAW-method J2Q]. During the opti- 
mization of the lattice constant we applied twice as high plane 
wave cut-off and 12 3 Monkhorst-Pack K-point set H27I1 . For 
the 5 12-atom supercell we used the T-point that provided con- 
vergent charge density. We plugged the optimized geome- 
try into the CPPAW supercell plane wave code with PAW- 
method that provides the hyperfine tensors B28ll . We applied 



2 



[Hi] 




[111] d) 



l±l> 
10) 



B=0 B>0 




e) 



l~U> 



i-U>^ Q Mo,t> 



Conduction band 





-ff -t- 






Valence band 





a. 



♦ 10) 




I0,t> 10,1 > 



FIG. 1: (Color online) a) The structure of nitrogen- vacancy (N-V) 
defect in our particular working frame. The place of vacancy is de- 
noted by an empty circle, b) The schematic single particle picture 
of the m s — 1 high spin states in ground state ( 3 A2,gs) and ex- 
cited state ( 3 _B,es). c) The fine structure of the 3 A2 and S E states 
at room temperature due to spin-spin interaction. Zero-field split- 
tings are D gs =2.88 GHz (H, D gs =1.42 GHz (H. During optical 
excitations the fluorescence is predominantly active for the m 3 — 
ground state (j0)) due to a non-radiative intersystem crossing of the 
± 1) es states with the many-body singlet states (not shown here), 
d) Splittings of the es substates in the presence of the magnetic field 
(B). LAC is expected between |0) and | — 1) states, e) Simplified 
energy-level diagramm with including the hyperfine structure asso- 
ciated with 15 N nuclear spin states | f ) and | J.) in the case of LAC 
regime of the applied B-field. At LAC, precession frequency be- 
tween excited-state sublevels |0, J.) and | — 1, j) can lead to nuclear- 
spin flip, which can be transfered to the ground state through nonra- 
diative intersystem crossing (curved arrow). 



the same basis set and projectors in both codes yielding vir- 
tually equivalent spin density of the defects. Other techni- 
cal details are given in Ref. Il9j. The charge density distri- 
bution was analyzed by the Bader-method H29fl . We briefly 
mention here that we provide the principal values of the hy- 
perfine tensors, called, An, A22, and A33 that can be found 
by diagonalization of the hyperfine tensors. If the hyper- 
fine field has C^ v symmetry then An=A22=A± and A33=Am 
where 1 1 means that hyperfine field coincides with the symme- 
try axis of the defect. The Fermi-contact term (a) is defined 
as a = [Aii + 2A±)/3 while the dipole-dipole term (6) as 
b = (Au — A±)/3. The hyperfine field is isotropic when 
6 = 0, i.e., Ax=A\\. 

First, we discuss the geometry of the defect. The obtained 
distance between the carbon atoms is 1 .54 A in perfect dia- 
mond. In the nitrogen-vacancy defect there are three carbon 
atoms (C a ) and one nitrogen atom (N) closest to the vacant 
site each possessing a dangling bond (Fig. QJ). The defect 
conserves its C3,, symmetry during the outward relaxation, 
and the dangling bonds of the C a and N atoms will point to 



the vacant site. The symmetry axis will go through the va- 
cant site and the N-atom which is the [111] direction in our 
particular working frame (see Fig.QJ). We found in our PBE 
calculations that the C a atoms are closer to the vacant site 
(1.64 A) than the N-atom (1.69 A) in the ground state. We 
obtained 5.97e Bader-charge on the N-atom that it is 0.97e 
larger than the number of valence electrons of the neutral N- 
atom. N-atom is more electro-negative than the C-atom. In- 
deed, the three C-atoms bound to N-atom (Cb) have 3.69e 
Bader-charge, so there is significant charge transfer from C& 
atoms toward the N-atom (~0.93e as total). That is the main 
source of the negative polarization of the N-atom. We found 
that the negative charge is distributed on many atoms around 
the defect. The C a atoms are even slightly positively polar- 
ized (3.97e) that will finally induce a dipole moment in the 
defect. Next, we briefly discuss the spin density of the ground 
state. The spin density is primarily originated from the un- 
paired electrons on the e defect level in the gap. Due to sym- 

J I I L 

metry reasons [ 19, 30] the e-level is only localized on the C a 
atoms but not on the N-atom. Therefore, large spin density 
is expected on the C a atoms while negligible on the N-atom. 
Indeed, very small hyperfine splitting was found for 14 N lEi 
and 15 N jl5lll8ll . However, the sign of the Fermi-contact term 
of the hyperfine interaction for 15 N was contradictory. Rabeau 
et al. assumed a negative value 11181 1 while Felton et al. has 



recently proposed a positive value [15j. The gyromagnetic 
factor of 15 N (7n) is negative, thus negative(positive) Fermi- 
contact hyperfine value (a) indicates positive(negative) spin 
density on the N-atom (n s ) because a ~ Tn^s- In our pre- 
vious LDA calculation 11911 we detected negative spin density 
on N-atom. Our improved PBE calculation justifies this sce- 
nario. Due to symmetry reasons the direct spin polarization 
of N-atom does not occur in the ground state but the large 
spin density on the C a atoms can polarize the core electrons 
of the N-atom, i.e., it is an indirect and weak spin polarization. 
In our PBE calculation we can also detect a slight dipole- 
dipole interaction for 15 N which is outside of our error bar 
(~0.3 MHz [ 19]). Our conclusions agree with the findings of 
Felton et al. [1 5] : the 15 N has positive hyperfine splittings and 
it is slightly anisotropic (see TableQ]i. The calculated hyperfine 
tensors for C a atoms agree nicely with the recent experimental 
values recorded at low temperature lfl5ll . The fair qualitative 
agreement between the experiment and theory for the ground 
state allows us to study the less known excited state with our 
tools. 

In the excited state we found a significant re-arrangement 
of atoms compared to the ground state. The C a atoms now 
farther from the vacant site (1.70 A) than the N-atom (1.63 A) 
in contrast to the case of the ground state. The N-atom attracts 
less electrons from its neighbor C-atoms: the Bader-charge of 
N-atom (5.88e) is 0.09e less than in the ground state. Conse- 
quently, the Bader-charge of Cb atoms increases by 0.03e. It is 
important to note that the Bader-charge of C Q atoms increases 
from 3.97e to 4.04e. Thus, the excitation induces a change in 
the dipole moment of the defect. Next, we discuss the change 
in the spin density distribution upon excitation. Fig. [2] shows 
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FIG. 2: (Color online) The calculated spin density difference be- 
tween the excited 3 E and the ground 3 A2 states of the NV cen- 
ter. We choosed two representative isosurface values indicated in 
the colored rectangles. The blue(cyan) ball(s) represent the nitro- 
gen(carbon) atoms. The vacant site is below the nitrogen atom. 



TABLE I: 3 A2 ground state (rows 2,3) versus 3 E excited state (rows 
4,5). The calculated principal values of the hyperfine tensor (columns 
2 to 4) compared to the known experimental data (columns 5 to 7) in 
MHz. The experimental data on 15 N is taken from Refs.QjlQl]. See 
text for the meaning of the question mark. 
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the calculated difference of the spin densities in the excited 
and ground states. As one can see that the spin density en- 
hanced a lot around N-atom (indicating with red lobes) while 
it dropped around the C a atoms (indicating with blue lobes). 
This can be explained by the hole left on the a\ defect level 
in the gap after excitation. The a\ defect level is significantly 
localized on the N-atom lfl9Tl . Thus, the spin polarization of 
the a\ defect level will spin polarize the N-atom considerably. 
Consequently, the spin polarization of the C a atoms will be 
smaller. According to the calculations the hyperfine constants 
of 13 C isotopes are dropped by around 50% (see Table J). 
However, the overall magnetization density of the N and C a 
atoms is about 95% the same both in the ground and excited 
states. In other words, the spin density mostly redistributed 
between the N-atom and the three C Q -atoms upon excitation. 
Indeed, it has been very recently found by using excited state 
spectroscopy in 15 N enriched diamond samples that the 15 N 
hyperfine signal is «20 times larger (~60 MHz) in the ex- 
cited state than in the ground state [ 16]. Our calculations can 
explain this feature. However, the applied model Hamilto- 
nian for describing the EPR of 15 N hyperfine signal (A (N) ) 
was incomplete in Refs. [l3. 16 . Fuchs et al. and Jacques et 



al. studied individual NV centers by confocal photolumines- 
cence miscroscopy where the actual defect was aligned to the 
[111] axis and the applied magnetic field was parallel to this 



axis, thus the angular dependence of A (N) was not measured. 
They assumed isotropic A (N) for the excited state while our 
study shows that it is strongly anisotropic. Fuchs et al. also 
noticed that A (N) ) should have the opposite sign in the ground 
and excited states 11611 - Our study shows that is posi- 
tive(negative) in the ground(excited) states in contrast to the 
previous assumptions Il3[|l6ll . 

Now, we discuss the consequence of our findings in the 
light of recent experiments on the dynamic polarization of 
single nuclear spins of the NV center BO. fl6ll . It has been 
demonstrated that the effective nuclear-spin temperature cor- 
responds to a fiK in this process [13] decoupled from the 
ambient temperature that can be the basic physical process 
in the measurement of zero-point fluctuations |2|]. In these 
measurements the de -polarization of the nuclear spins of 15 N 
El 1611 and 13 C a II 1311 have been demonstrated. This has been 
achieved by the level anticrossing (LAC) of the electron spin 
m s sublevels in the excited state. The LAC effect may ap- 
pear if the m s sublevels cross at a given external magnetic 
field (see Fig.[T];,d). We show a refined and corrected model 
of Ref. 13 accounted for LAC. We study the de-polarization 
of 15 N isotope but it can be generalized for the 13 C isotopes 
straightforwardly. The Hamiltonian of the system (with ne- 
glecting the nuclear-Zeeman splitting) can be written as Jl3ll . 



H = D es S 2 z + g ef i B BS z 



A es SI , 



(1) 



where S_ and J are the electron and nuclear-spin operators, 
D es the excited-state zero-field splitting, g e the electron g 
factor, 11 b the Bohr-magneton, and A es the excited hyper- 
fine coupling. We assume positive B-field. Because A es is 
anisotropic, the A es S_ J term can be written with the a es and 
bes hyperfine splittings and the spin-shift operators as, 



S+ 1- + S- I 4 



V^es ^es 

) + S z I z (a es + 2b es ) (2) 



The hyperfine field of 5 N is parallel to the symmetry axis, and 
(a es - b es )=A±fa-39 MHz while (a es + 26 es )=A||W-58 MHz. 
According to a recent study llzjl D es =+1.42 GHz (m s = 
sublevel is below m s = ±1 sublevels), so we can restrict our 
study to the excited state m s = and m s = — 1 sublevels 
(see Fig. HJl). In the basis [| - 1,4); | - 1, f); |0, |); |0, t); ] 
and by choosing the origin of energy level at level |0, \), the 
Hamiltonian described by Eqs. lll2l can be written as 



H = 
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c d 
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7J es ±A||/2 
c = g e [i B B 
d = A ± /V2 



The eigenstates of this Hamiltonian are |0, |), | — 1,4), 
|+) = a|0, 4) + 0\ - 1, T) and |-> = 3\0, 4) - a\ - 1, f). 



By following the arguments in Ref. 13, the transition from 
the ground state |0, f) to the excited state remains nuclear 
spin conserving, whereas the transition from |0, 4) results in 
(a|+) + 0\— )) in the excited state (see Fig.OJ). This super- 
position state then starts to precess between the appropriate 
states at frequency Q = l/{2h) x - c) 2 + 4d 2 ] 1/2 , 
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where h is Planck's constant. The precession frequency de- 
pends on B via electron Zeeman-effect (c in our notation) that 
is minimal at LAC resonance, and will be equal to \d\/U = 
\A±\/ \/2h. Jacques et al. assumed isotropic hyperfine split- 
ting for 15 N, therefore, they applied «60 MHz in this formula 
lfl3[| . Our analyzis shows that rather the | A± | «39 MHz should 
be substituted here. Nevertheless, this precession frequency is 
still at the same order of magnitude as the excited state decay 
rate, 12 ns fl^l . Thus, the spin-flip process is very efficient 
between |0, J.) and | — 1,1) states, and we can exp lain the 
de-polarization of 15 N found in the experiments Gamin. 
It has been found that the probability of the de-polarization 
effect significantly depends on the misalignment of the mag- 
netic field from the symmetry axis 01311 . This may be partially 
explained by the anisotropy of the 15 N hyperfine splitting be- 
side the mixing of the spin states. 

We found other intriguing properties of the spin density in 
the excited state. As apparent in Fig. |2j the spin density and 
the corresponding hyperfine tensors change considerably also 
for the proximate 13 C isotopes. We show the hyperfine tensors 
only for the most significant change, for C a atoms in Table U 
Beside that new 13 C isotope becomes active above N-atom 
that has negligible spin density in the ground state. In addi- 
tion, the spin density of the sets of 6(3) C-atoms at R=3.9 A 
decreases(increases) due to excitation where R is the distance 
from the vacant site. According to our previous study M19I1 . 
one of these 13 C isotopes was manipulated in the qubit and 
quantum register applications yl [ljj] . Our study shows that 
during the optical set and read-out processes the spin-density 
of the addressed proximate 13 C isotopes changes indicating an 
effective oscillating magnetic field with the inverse lifetime of 
the excited-state. This may also influence the decoherence of 
the entangled electron-nuclear spin state that has not yet been 
considered 
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